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Abstract 

The sensitivity of optimal integrated observables to electroweak triple gauge 
couplings is investigated for the process e~e"'" — ^ W~W^ — >■ 4 fermions at fu- 
ture linear colliders. By a suitable reparameterisation of the couplings we 
achieve that all 28 coupling parameters have uncorrelated statistical errors 
and are naturally normalised for this process. Discrete symmetry properties 
simplify the analysis and allow checks on the stability of numerical results. 
We investigate the sensitivity to the couplings of the normalised event dis- 
tribution and the additional constraints that can be obtained from the total 
rate. Particular emphasis is put on the gain in sensitivity one can achieve 
with longitudinal beam polarisation. We also point out questions that may 
best be settled with transversely polarised beams. In particular we find that 
with purely longitudinal polarisation one linear combination of coupling pa- 
rameters is hardly measurable by means of the normalised event distribution. 
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1 Introduction 



The Standard Model (SM) of electroweak interactions has been thoroughly investi- 
gated both theoretically and experimentally and turned out to be very successful. 
After the discovery of the massive gauge bosons W and Z, the direct measurement 
of the triple gauge couplings (TGCs), i.e. the couplings between two charged and one 
neutral gauge boson, has been an important issue at the TEVATRON and LEP as 
their values are determined by the non-Abelian electroweak gauge group. At e~^e~ 
colliders the production of W pairs, single Ws, single photons and single Zs is suitable 
for that. In this paper we consider the process e~e~^ W~W^ /1/2/3/4, where 
both the jWlV and the ZWW couplings can be measured at the scale given by the 
cm. energy. 

Deviations from the SM at the 'yWW and ZWW vertices can be parameterised 
in a general framework. Allowing for complex couplings the most general vertex 
functions lead to altogether 28 real parameters . All four LEP collaborations have 
investigated TGCs 0, the tightest constraints being of order 0.05 for Agf and A^, 
of order 0.1 for Ak^ and of order 0.1 to 0.6 for the real and imaginary parts of 
C and/or P violating couplings. All these values correspond to single parameter 
fits. Since only up to three-parameter fits have been performed, only a small subset 
of couplings has been considered at a time, thereby neglecting correlations between 
most of them. Moreover many couplings, notably the imaginary parts of C and P 
conserving couplings, have been excluded from the analyses. 

At a future linear e~^e~ collider like TESLA or CLIC @], respectively covering 
a cm. energy range from about 90 GeV to 800 GeV and from about 500 GeV to 
5 TeV, one will be able to study these couplings with unprecedented accuracy. In 
this way one may begin to be sensitive to different extensions of the SM, which 
predict deviations at the TGCs from the SM values, typically through the effects of 
new particles and couplings in radiative corrections. Some examples, where effects of 
order 10~^ may occur, are supersymmetric models models containing several 

Higgs doublets 0, §], Eq vector leptons ^ or Majorana neutrinos [|10| and the minimal 
3-3-1 model [|1T|. For left-right symmetric models [^2|, |13|, ^ and mirror models |13 



the effects are predicted to be much smaller, whereas models containing composite 
W bosons or an additional gauge boson Z' |T3| may lead to larger effects. 

Given the intricacies of a multi-dimensional parameter space, the full covariance 
matrix for the errors on the couplings should best be studied. The high statistics 
needed for this will for instance be available at TESLA, where the integrated lu- 
minosity is projected to be about 500 fb~ or more per year at 500 GeV which for 
unpolarised beams amounts to about 3.7 million produced W pairs (without cuts). 
For a run at 800 GeV the luminosity is expected to be twice as high, leading to 
3.9 million W pairs. Moreover, polarised beams will be particularly useful to disen- 
tangle different couplings. In fact, certain directions in the parameter space of the 
couplings, the so called right handed couplings, are difficult to measure in W pair 
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production with unpolarised beams II^]. In this case their effects are masked by the 
large contribution from the neutrino exchange. With polarised beams the strength of 
the neutrino exchange contribution can in essence be varied freely. 

In experimental analyses of TGCs and various other processes, optimal observ- 
ables 16 1 have shown to be a useful tool to extract physics parameters from the 
event distributions. These observables are constructed to have the smallest possible 
statistical errors. Due to this property they are also a convenient means to deter- 
mine the theoretically achievable sensitivity in a given process. In addition, they take 
advantage of the discrete symmetries of the cross section. 

The method proposed in |18] allows one to simultaneously diagonalise the covari- 
ance matrix of the observables and the part of the cross section which is quadratic 
in the couplings. In this way one obtains a set of coupling constants that are nat- 
urally normalised for the particular process and — in the limit of small anomalous 
couplings — can be measured without statistical correlations. This allows one to see 
for which directions in parameter space the sensitivity to the TGCs is high and for 
which it is not, which is hardly possible by looking at covariance matrices of large 
dimension without diagonalisation. Moreover, the total cross section acquires a par- 
ticularly simple form and provides additional constraints. 

The basic purpose of this work is to use this extended optimal observable method 
for a systematic investigation of the prospects to measure the full set of TGCs in 
W pair production at linear collider energies, with special emphasis on initial-state 
polarisation. The usefulness of the method becomes particularly evident when consid- 
ering imaginary parts of CP conserving TGCs. In this subspace of couplings we find 
one direction to which — in the linear approximation and for longitudinally polarised 
beams — there is no sensitivity in the process we consider. In this was not taken 
into account and led to numerical instabilities. It is therefore essential to disentangle 
the measurable TGCs from the hardly measurable one. As a historical motivation one 
may think of the electromagnetic nucleon form factors Fi and F2, where the choice of 
linear combinations Ge and Gm leads to a simplification of the Rosenbluth formula 
for the differential cross section of electron- nucleon scattering (see e.g. |]19[)- Since 
we deal with 28 couplings here, an appropriate choice of parameters is essential. 

We restrict ourselves to the semileptonic decays of the W pair, where one W bo- 
son decays into a quark-antiquark pair and the other decays leptonically, but leave 
aside the decay into tUt- since these events have a completely different experimental 
signature. The selected channels have a branching ratio of altogether 8/27, which is 
six times larger than that of both Ws decaying into ez/e or fiu^. On the other hand. 



the channels where both Ws decay hadronically are difficult to reconstruct [g^]- The 
semileptonic channels have only one ambiguity in the kinematical reconstruction if 
the charges of the two jets from the hadronically decaying W are not tagged. Then 
one cannot associate the jets to the up- and down-type (anti)quark of the W decay, 
and therefore has access only to the sum of the distributions corresponding to the 
two final states. 
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This work is organised as follows: In Sect. we recall the helicity amplitudes and 
cross sections of the process using a spin density matrix formalism. In Sect. |^ the 
optimal-observable method is presented in the form as it is used in our numerical 
calculations. We explain our technique to implement the simultaneous diagonalisa- 
tion in a numerically stable way. The role of discrete symmetries in the framework 
of optimal observables is described. Many of the symmetry relations are well known, 
notably the classification of the TGCs into four symmetry classes and its appli- 
cability to the optimal-observable method [|16|, |1^]. Other properties are used for a 
check on the numerics. In Sect. ^ the dependence of the sensitivity on longitudinal 
beam polarisation is illustrated by a simple model. In Sect. |^ we show analytically 
that one is insensitive to one of the imaginary CP conserving couplings in the case of 
longitudinal beam polarisation. However, this particular coupling becomes accessible 
with transverse beam polarisation. In Sect. |^ we present our numerical results, in 
Sect. 0our conclusions. 

2 Cross Section 

First we briefly recall the differential cross section of the process 

/i + /2 ^h + f, 

for arbitrary initial beam polarisations, where the final state fermions are leptons or 
quarks. Our notation for particle momenta and helicities is shown in Fig. 0. Our 
coordinate axes are chosen such that the e~ momentum points in the positive z- 
direction and the y unit vector is given by = (k x q)/|k x q|. 




Figure 1: Momenta and helicities of the particles in the e'^e cm. frame. 
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In the e^e cm. frame and at a given cm. energy a/s, a pure initial state of 
longitudinally polarised e~ and e"*" is uniquely specified by the e~ and e'^ helicities: 

|rr) = |e-(k,T)e+(k,T)) (t,t = ±1). (2) 

Note that fermion helicity indices are normalised to 1 throughout this work. A mixed 
initial state of arbitrary polarisation is given by a spin density matrix 

p=Y, \'rr)p{rr){T'r'){T'T\, (3) 
(r) 

where (r) denotes summation over r, r, r' and r', and the matrix elements satisfy 
P(^:5r^(^/:5r/-) = P(t't'){tt) ^ud J2t,t P(tt)(t-) = 1- Wc dcfiuc the cross section operator 

da- ^^\r'r')da^r'T')irT){rr\ (4) 

by requiring the differential cross section for arbitrary p to be 

da\p = tr(do-p) = ^ (i(T(r/:f/)(^:f) (5) 

The matrix da(^r'T'){TT) is given by 

C^^(r'r')(rr) = ^ / (/|T|rr) (/ITItV)* , (6) 

where we neglect the electron mass in the flux factor. Here T is the transition 
operator, |/) = |/i (pi, ti) /a (p2, T2) /s (ps, T3) /4 (p4, T4)) the final state and 

the usual phase space measure for final states. Using the narrow-width approximation 
for the Ws the result is 

R ^ ^ fg) 



Here m^^ is the W boson mass, Vw its width and {3 = ^1 — irri^r/s its velocity in 
the e'^e~ cm. frame. The sum (A) runs over A, A', A and A'. Furthermore 



dvl^^j,\^^^^ = d(cose) (i$ (AA, e\T\TT) (A'A', e|r|r¥)* (10) 
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Figure 2: Definition of azimutfial angles. 



is tlie differential production tensor for the W pair and 

dVyx = d{cos^)d^{hj2\r\x){fj;\r\\')\ 

dVyj = d{cos^)dip {fjl\r\\) {hh\T\\r (11) 

are the tensors of the W~ and decays in their respective cm. frames. Note that 
in the narrow-width approximation the intermediate are treated as on-shell. We 
define the W helicity states which occur on the right hand side of (p!OD in the frame of 
Fig. |T|, i.e. we choose the e~e^ W~W^ scattering plane as the x-z-plane and define 
to be the polar angle between the W~ and e~ momenta. We choose a fixed direc- 
tion transverse to the beams in the laboratory. By $ we denote the azimuthal angle 
between this fixed direction and the e~e+ — > W~W^ scattering plane (see Fig. 0(a)). 
The respective frames for the decay tensors (|lID are defined by a rotation by about 
the y-axis of the frame in Fig. |I| (such that the W~ momentum points in the pos- 
itive z-direction) and a subsequent rotation-free boost into the cm. system of the 
corresponding W . The spherical coordinates are those of the /i and momentum 
directions, respectively. In its rest frame, the quantum state oi a.W boson is deter- 
mined by its polarisation. For real Ws we take as basis the eigenstates of the spin 
operator Sz with the three eigenvalues A = ±1,0. For off-shell Ws a fourth, scalar 
polarisation occurs but is suppressed by mf/mw in the decay amplitude. 
Neglecting the electron mass we have in the SM 

dai^r'-')(T-) = for r = r or r' = t', (12) 

which we will use in the following (this point is further discussed in Sect. p.3| ). At 
a linear collider the initial beams are uncorrelated so that their spin density matrix 
factorises, i.e. 

P(t-){t'-') = Ptt'PTt' , (13) 

where prr' and 'p^:^, are the two Hermitian and normalised spin density matrices of 
e~ and e"*" respectively. We parameterise these matrices by 
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Figure 3: Feynman diagrams for the process e with anomalous TGCs. 
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where < y?^ < 27r, and the vector components of a are the Pauh matrices. The 
degrees of transverse and P(^ of longitudinal polarisation obey the relations 
{P^f + {P^f < 1 and Pt > 0. The components of in (|TB|) refer to the frame of 
Fig. |l]. Note that choosing the same reference frame for p~ and p~^ while specifying 
each spinor in its respective helicity basis results in different forms of the density 
matrices in (|1^. The relative azimuthal angle ip = (f~ — ip^ between p~ and p^ is 
fixed by the experimental conditions, whereas the azimuthal angle ip~ of p^ in the 
system of Fig. |1| depends on the final state (see Fig. 0(b)). For the case where P^" ^ 
we choose the transverse part of the vector p~ as fixed direction in the laboratory. 
Then we have $ = — v?~- Using (^ and (|T2]) to (|T3|), we obtain the differential cross 
section for arbitrary polarisation: 



dcT\, = J|(i + pr)(i-^^^)^^(+-){+-) 

+ (l-Pr)(l + P,+)rfa(_+)(_+) 



2P 



Pt 



Re (i(T(+_)(_+) cos {ip + 2$) 
-){-+) sin (^ + 2$) 



+ Im da I 



{+- 



(16) 



In the absence of transverse polarisation, (|T6D is independent of $ due to rotational 
invariance. 

In our analysis we evaluate the matrix elements in (^) at tree level in the SM, re- 
placing the 'jWW and ZWW vertices by their most general forms allowed by Lorentz 
invariance. The corresponding Feynman diagrams are shown in Fig. ^. New physics 
may also lead to deviations from the SM values at the fermion-boson vertices . 
For these vertices we however retain the SM expressions, following the rationale that 
at present they are confirmed experimentally to a much higher precision than the 
triple gauge couplings. We remark that there are scenarios of physics beyond the 



8 



SM where such a treatment is not adequate, since the process e^e+ — > W~W~^ can 
receive non-standard contributions that cannot be expressed in terms of anomalous 
fermion-boson or three-boson couphngs (an example are box graph contributions in 
super symmetric theories ^). Such effects can still be parameterised within a 



more general form factor ansatz |23]. If they are important, an analysis in terms of 
only TGCs will not give a correct picture of the underlying physics, but it will still 
correctly signal a deviation from the SM. We finally remark that radiative correc- 
tions in the SM itself can be included in the analysis procedure we develop here (see 
Sect. ^). The purpose of the present study is however to investigate the pattern of 
sensitivity to TGCs and its dependence on beam polarisation, and for this purpose 
it should be sufficient to take the SM prediction at tree level. 

For the three-boson vertices we use the parameterisation (2.4) in Ref. [|I|, which 
we express in terms of the coupling parameters , Ky, Ay, , g^ , and Xy 
{V = 7, Z) by means of the transformation (2.5) in the same reference. In other words, 
we understand these parameters as form factors of the three-boson vertices, which 
depend on the boson virtualities and can take complex values, and not as coupling 
constants in an effective Lagrangian, which by definition are energy independent and 
real- valued. 

For a given e~ beam helicity r the process (|l]) is not sensitive to all couplings, but 
only to the linear combinations gi, hl, etc. for left handed (r = —1) and g^, kh, etc. 
for right handed (r = 1) electrons, where 

g^ = Asm' ewgj + {2 -4sm^ew)^gf, 

g^ = Asm'ewgJ-Asin'ew^gf, (17) 

and similarly for the other couplings |l|, [T^. Here e denotes the positron charge, 6w 
the weak mixing angle, and ^ = s/ (s — m|) the ratio of the Z and photon propagators. 
The parameterisation (^) will in the following be called the L(eft)-R(ight)-basis. 

The expressions of the amplitudes can be found in |l|]. For convenience of the 
reader we rewrite the W~W^ production part in terms of the LR-basis. The matrix 
element of (|10D is given by 

(AA,e|T|rr) = -v^e^>[(r; A, A; 9) r?dt_^,(e), (18) 

where r/ = Ar(-l)^, Ar = |(r - r), AA = A - A, and Jo = max(|Ar|, |AA|). The 
definition of the d-functions and our spinor conventions, as well as the SM matrix ele- 
ments for the W decays in (|ll|) are listed in the Appendix. The production amplitude 
M. = Ai^'^'^ + Ai'^ consists of two terms given by 
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Table 1: Coefficients A'^^j, B^j and C-^j of the matrix elements (|TP|) and (PP]). The 
indices of the couplings are a = L for r = — 1 and a = R for r = +1. The relativistic 
factors are defined by 7 = ^/s/2mw and /5 = (1 — 4m^/s)^/^. 



(AA) 






^AA 


(+-),(-+ 


) 





2^2/3 


(+0) 
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27 
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(++) 
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1 


1/7^ 


(— ) 


gl + 272A, - ip-^ka + + 27'/5)A. 


1 


1/7^ 


(00) 


^? + 27'fi:a 


27^ 


2/7^ 



The expressions for (which contains the left handed couplings for r = — 1 and the 
right handed ones for r = +1) and for Byj;^ and C;^^ are listed in Table |I]. Since the 
vector bosons carry spin 1, the TGCs do not contribute to the WW helicity combi- 
nations (H — ) and ( — h). For the other helicity amplitudes, the largest power of 7 in 
the coefficients AZ-^ coincides with the number of longitudinal Ws. An exception are 
the couplings Aa and Aa, which correspond to dimension-six operators in the effective 
Lagrangian (cf. and occur with an additional factor of 7^. Note that the largest 
kinematical factors in A\-^ behave like 7^ at high energies, in contrast to the basis of 

AA 

form factors used in Table 4 of where huge factors of 7 appear. In the SM at 
tree level one has 

gX = l, Kv = l {y = l.Z) (21) 

and all other couplings equal to zero. For the anomalous parts of the couplings we 
write !\g( = g( — \ and A/ty = — 1 as usual. 

A detailed discussion of the differential cross section is given in [|l|. Here we 
only point out some salient features of the high-energy limit. In Fig. ^ we show the 
total cross section for unpolarised beams as a function of y^. It rises rapidly from 
threshold up to a maximum of about 20 pb at a/s ~ 200 GeV, and in the SM decreases 
for higher cm. energies. In the SM each Z-coupling is equal to the corresponding 
photon coupling. Since ^ = 1 + 0(7"^) the L-couplings are then of order 1 and the 
R-couplings of order 7"^. For r = +1 this leads to a high energy behaviour of at 
most M. ~ 0(1). For r = — 1 we have (^f = = 2^ in the SM, and the coefficients 
Ayj^ and -B^^j only differ by a factor 2 + 0(7"^) according to Table |l|. As they occur 
with different sign in (^) and ( pUj ) this again results in a high-energy behaviour 
M. ~ 0(1), except for very forward W~ momentum where there is an enhancement 
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Figure 4: Total cross section with unpolarised beams as a function of the cm. energy 
in the SM and for one anomalous coupling differing from zero. Some curves coincide 
as explained in the text. 

by the propagator factor (1 + — 2/3 cos 6)~^. Altogether, these gauge cancellations 
preserve the unitarity of the SM. We also plot in Fig. ^ the total cross section for 
one anomalous coupling differing from zero. At high energies each coupling mainly 
contributes via the W helicity amplitude where it occurs with the highest power 
of 7, i.e. either linearly or quadratically according to Table |l|. At sufficiently high 
energy, the square of an anomalous term dominates over its interference term with 
the SM amplitude. In the limit (3 —>■ 1 the couplings gi, g^, g^ and k enter with 
a factor 7, whereas n, A, A enter with a factor 7^, which explains their different 
behaviour in Fig. ^ Some couplings have equal coefficients in this limit, which leads 
to a degeneracy of the curves. We also remark that even if more than one anomalous 
coupling differs from zero, anomalous amplitudes belonging to couplings of different 
C or P eigenvalue do not interfere in the total cross section with unpolarised beams 
(cf. Sect. U). 



^ , , , I , , , , I , , , , I , , , , I , , , , 
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- AgJ, gj, gl or k-^ = 0.1 

' Agf , gf , gl or Kz = 0.1 ' 
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r = 0.1 ... - 

I A.y or A., = 0.1 I 

- or = 0.1 ..--''1.-'' ^.---"'"^ - 
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3 Optimal Observables 



In this section we discuss how the method of optimal integrated observables is applied 
to our case. In an experiment one measures the differential cross section 

S = da\p/d(j), (22) 

where (p denotes the set of all measured phase space variables. We distinguish be- 
tween the information from the total cross section a = J da and from the normalised 
distribution S/a of the events. We first investigate how well TGCs can be extracted 
from the latter, and then use a to get constraints on those directions in the space of 
couplings to which the normalised distribution is not sensitive. 

Since 5 is a polynomial of second order in the anomalous couplings (defined as 
the couplings minus their tree-level values in the SM) one can write 

S = Sq + Y^ Siihi + ^ S2ijhihj . (23) 

Here 5*0 is the tree-level cross section in the SM, and each hi is the real or imaginary 
part of an anomalous TGC Thus there are 28 real parameters hi altogether. In the 
remainder of this section i and j each run from 1 to 28. 

One way to extract the coupling constants from the measured distribution (|23|) is 
to look for a suitable set of observables Oi{(f)) whose expectation values 

E[0,] = - fdcpSO, (24) 



cr 

are sensitive to the dependence of S on the couplings hi. Given the present exper- 
imental constraints, we expect these variations to be small. Thus we expand the 
expectation values to first order in hf. 

E[Oi] = Eo[0.i] + J2ci,h, + 0{h''), (25) 

j 

with 

Eo[0,] = — fd<j)SoO, , (26) 

(To 

c = L fdcf) OiSi, [d(P SoOi , (27) 

= Jd(j) Sij . (28) 

Here -EolC^t] is the expectation value for zero anomalous couplings, whereas Cij gives 
the sensitivity of E[Oi] to hj. Solving ( P^D for the set of the hj we get estimators for 
the anomalous couplings, whose covariance matrix is given by 

V{h) = j^c-'V{0){c-Y, (29) 
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where we use matrix notation. Here is the number of events, and 
V{0)^^ = - fd^SoO,0,-Eo[0,]Eo[0,] + 0{h) 

■' (Tn J 



(30) 



is the covariance matrix of the observables, which we have Taylor expanded around 
its value in the SM. As observables we choose 



Su{4>\ 

Som^ 



From (|27|) and ( pOD one obtains for this specific choice 

V{0)=c + 0{h), 



and therefore 



0{h). 



(31) 



(32) 



(33) 



The observables (0) are "optimal" in the sense that for hi ^ the errors ( p3D on 
the couplings are as small as they can be for a given probability distribution.^ Apart 
from being useful for actual experimental analyses, the observables ( ^I]) thus provide 
insight into the sensitivity that is at best attainable by any method, given a certain 
process and specified experimental conditions. In the case of one parameter this 
type of observable was first proposed in |]T^, the generalisation to several parameters 
was made in |T^. Moreover, it has been shown that optimal observables are unique 
up to a linear reparameterisation [Q. We further note that phase space cuts, as 
well as detector efficiency and acceptance have no influence on the observables being 
"optimal" in the above sense, since their effects drop out in the ratio (|3l|) . This is 
not the case for detector resolution effects, but the observables (^1]) are still close 
to optimal if such effects do not significantly distort the differential distributions Su 
and 5*0 (or tend to cancel in their ratio). To the extent that they are taken into account 
in the data analysis, none of these experimental effects will bias the estimators. 

In the present work we use the method of optimal observables in the linear approx- 
imation valid for small anomalous couplings. But we emphasise that the method has 
been extended to the fully non-linear case where one makes no a priori assumptions 
on the size of anomalous couplings in [O]. Such a non-linear analysis for real data 



was presented in and turned out to be very convenient and highly efficient. 

Given the projected accuracy at linear colliders, it will in general be necessary 
to take into account radiative corrections to the process e~e+ — > /1/2/3/4 within the 

p^] . One possibility to 
"deconvolute" these 



SM, which have been worked out in detail in the literature 
include them in searches for non-standard TGCs would be to 
corrections. For this write the one-loop corrected differential cross section in the SM 
as 

SsM,coA<P)= /rf0'5o(0')F(0',0), (34) 



*For details on this so called Rao-Cramer-Frechet bound see e.g. 
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where 5*0 is the tree-level expression, and the integral kernel can e.g. be obtained 
from an event generator by generating events according to both S'sM.corr and 5*0. 
Approximating the true physical cross section as 

Vys(0)= fd^'S{4>')F{(P',(P), (35) 



with S given as in (23), one could invert this convolution bin by bin, and then extract 



the couplings hi from the deconvoluted Born level cross section S as described before. 
The error made in (^) is that the SM radiative corrections encoded in F do of course 
not apply to the anomalous part 5" — 5*0 of the cross section, but this error is of order 
hi times the weak coupling constant ctw Should effects beyond the SM be found 
in such an analysis, one would in a second step have to consider more sophisticated 
methods to quantitatively disentangle them from SM radiative corrections. 

3.1 Simultaneous Diagonalisation 

Although discrete symmetry properties (see Sect. ^73|) reduce the 28x28 matrices 



V{h), c, etc. to blocks of 8x8 and 6x6 matrices, these blocks still contain many 
non-negligible off-diagonal entries. To make explicit how sensitive the process is to 
each direction in the space of couplings and to identify the role of polarisation we 
need to know directions and lengths of the principal axes of the error ellipsoid defined 
by (ppp, i.e. we have to know the eigenvalues and eigenvectors of V{h). Using optimal 
observables, to leading order in the hi, the three matrices V{h), c and V{0) are 
automatically diagonalised simultaneously due to (|32D and (^3|) . This means that 
after such a transformation each observable is sensitive to exactly one coupling, and 
the observables as well as the estimators of the couplings are statistically independent. 
Since V{h) is symmetric the diagonalisation could be achieved by an orthogonal 



transformation. Following the proposal of [|18[ we take however a different choice and 



transform simultaneously V{h) into diagonal form and the normalised second-order 
part of the total cross section into the unit matrix: 

(^2ij = — fdcj) S2ij Sij . (36) 

(To -J 

This can always be done since d'2 is symmetric and positive definite. We therefore 
arrive at the following prescription for the transformation of the couplings (using 
vector and matrix notation): 

h h' = A^h, (37) 
V(h)-' Vih')-' = A^V{h)-'A = diag {{5h[)-^ {5h',)-\ . . . , {5h',,r') , (38) 

(J2 ^2 = A^o-gA = 1, (39) 

where 5h'i are the one-sigma errors on the new couplings. This transformation exists 
and is unique up to permutations and a sign ambiguity for each h[. Note that the 
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matrix A is in general not orthogonal. From (p3D, (|3T|) and (|37|) the transformation 
of all other quantities follows as 

^1 ^ 5; = A^5i, 

O 0' = A^0, 

c c' = A^cA, 
V{0) V{0') = A^V{0)A. (40) 

The meaning of (|39D is that all quadratic terms contribute to the total cross section 
with equal strength: 



\ i=i 1=1 / 

where a'^ = a^^ jdcfyS'^. Thus the anomalous couplings do not mix in a and are 
"naturally" normalised for the process which we consider. This is not true in the 
conventional basis, where changing different anomalous couplings by the same amount 
has completely different effects on the total cross section (see Fig. ^ in Sect. ||). 

Moreover, the particularly simple form (^) of a easily allows one to use the infor- 
mation from the total rate: it constrains the couplings to lie between two hyperspheres 
in the space of the h'^, whose difference in radius depends on the measurement error 
on a. Making in addition use of the sign ambiguity in ([37|), one can for all ^ 
choose the sign of h[ such that a'^ > 0. This choice is however not relevant for the 
analysis. 

We finally note that the presented method of simultaneous diagonalisation is quite 
similar to the way one analyses the normal modes of a mult i- dimensional harmonic 
oscillator in classical mechanics p9| . There the harmonic potential (corresponding to 



V) is diagonalised with respect to the scalar product defined by the kinetic energy 
(corresponding to (T2). 

3.2 Numerical Realisation 

We now give some details of how the simultaneous diagonalisation can be carried 
out numerically. Although the procedure finally aims at the disentanglement of the 
couplings as achieved in (|38D, the numerical computation of V{h) or V{h)~^ from ( pQ] ) 
needs the inverse of the matrix c or V{0) and might therefore be unstable, because 
before the diagonalisation one cannot single out those directions in parameter space 
where the errors are large. However, as V{h') and V{0') are simultaneously diagonal 
for our observables, we can according to ( |5^ ) equally well compute the diagonal entries 
of 

V-(0')=c' = diag(c;,4,...,48), (42) 
and extract the errors on the couplings using 



6K = (43) 
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Hence, using the shorthand notation V = V{0) and V = V{0'), we have to solve 
the + n equations 



A^asA = 1, 

A^VA = V (44) 
for the 71^ entries of A and the n diagonal elements of V . Since the multiplication of 



S* by a constant changes neither the observables (0) nor V nor o"2, the matrices A 
and V only depend on the normalised distribution of the events but not on the total 
rate N . The latter enters the errors on the transformed couplings only through the 
statistical factor N~^^'^ in (|^). From (|^) it follows that 

VA = a2A v. (45) 

This is a generalised eigenvalue problem, with the c[ being the generalised eigenvalues 
and the columns of A being the generalised eigenvectors. The pair (V , A) is called 
the "eigensystem" of (^If). A standard method for solving (^) is to first perform a 
Cholesky decomposition p7| , pB[ 



<T2 = MM^, (46) 

where M is a lower triangular matrix, i.e. Mij = for j > i. Algorithms for the 
computation of M can be found in the same references. Then (^5)) is equivalent to 



CX = XV', (47) 

where 

C = M-W{M-^f, (48) 

X = M^A. (49) 

Equation (|^) denotes a (usual) eigenvalue problem for the matrix C, whose eigenval- 
ues are the same as the original ones and whose eigenvectors are the columns of X. 
Since C is symmetric, ( |T7| ) can be solved and the eigenvectors are orthogonal with 



respect to the standard scalar product (assuming non-degeneracy of the eigenvalues). 
Requiring the eigenvectors to be normalised to 1, we have n conditions 

X^X = 1, (50) 



which together with the equations (^Tf ) are equivalent to (@). The generalised 
eigenvectors are obtained by solving (|49|) for A. The procedure has to be followed for 
each initial state polarisation and cm. energy separately, leading in general to differ- 
ent eigenvalues and transformation matrices; the dependence of those quantities on 
polarisation is investigated in Sect. In each case we use the procedure iteratively. 
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i.e. once A is obtained, we compute V and a2 for the transformed observables and cou- 
plings, and then diagonahse these — aheady approximately diagonal — matrices again. 
We found this to be essential to assure the numerical stability of the results. A stable 
value was reached in most cases by the second evaluation and at the latest by the 
fourth. In all cases at least five evaluation steps were carried out. The numbers 
presented in Sect. ^ were obtained by averaging over the results of several subsequent 
steps where the stable value had already been reached. 

We note that for situations where o"2 and V have the same block diagonal structure, 
the diagonalisation can of course be carried out for each block separately. This is 
relevant in the presence of discrete symmetries. 



3.3 Discrete Symmetries 

Let us now discuss the special role of the combined symmetry operations CP and 
CPT in the context of our reaction |18|. Here C denotes charge conjugation, P 



parity reversal, and T "naive time reversal" , i.e. the reversal of all particle momenta 
and spins without the interchange of initial and final state. Under the condition that 
the initial state, as well as phase space cuts and detector acceptance are invariant 
under a CP transformation, a CP odd observable gets a nonzero expectation value 
only if CP is violated in the interaction. Similarly, if the initial state, phase space cuts 
and acceptance are invariant under CPT followed by a rotation by 180° around an 
axis perpendicular to the beam momenta, a nonzero expectation value of a CPT odd 
observable implies the interference between absorptive and nonabsorptive amplitudes 
in the cross section. In terms of the three-boson couplings one finds that to 0{h) 
the expectation values of CP even (odd) observables only involve the CP conserving 
(violating) couplings gi, k, A, {g4, k, A). Similarly, CPT even (odd) observables are 
to first order only sensitive to the real (imaginary) parts of the coupling parameters. 
The coefficient matrix c is thus block diagonal in the following groups of observables: 

(a) : CP and CPT everi, 

(b) : CP even and CPT odd, 

(c) : CP odd and CPT even, 

(d) : CP and CPT odd. 

One further finds that the first-order terms an in the integrated cross section can 
only be non-zero for couplings of class (a). 

The above requirements on the initial e^e~ state are nontrivial in the case of 
polarised beams. Since charge conjugation exchanges e"*" and e~, a CP invariant spin 
density matrix requires p~ = p^, in particular Pf = P^' and Pf = —Pi^- Under the 
same conditions the density matrix is also invariant under CPT times a rotation by 
180° around k x p~ . Let us investigate the situation where these requirements are not 
satisfied. From Fig. |]we see that the e~^e~ states where the beam helicities are aligned 
are CP eigenstates. But the antialigned (Ar = 0) states are interchanged under CP 
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Figure 5: Effect of a CP transformation on an e+e state witli fielicities aligned (top) 
or antialigned (bottom). 

and tlierefore are not eigenstates. On tlie otlier liand, tlie amplitudes for these helicity 
combinations are proportional to the electron mass in the SM (supplemented by the 
most general TGCs). They are thus generically suppressed by rrie/Mw compared 
with the amplitudes for aligned beam helicities (Ar = ±1).0 With transverse beam 
polarisation, the two types of amplitudes can interfere, giving small rrie/My/ effects in 
the cross section. With purely longitudinally polarised beams they do not interfere, so 
that effects due to the Ar = combinations of e"'"e~ are of order {rrie/MwY thus 
beyond experimental accuracy. We remark that the same holds e.g. in the minimal 
supersymmetric extension of the SM, where left- and right handed leptons as well as 
their superpartners only mix with a strength proportional to the lepton mass. To 
investigate what can happen in more generic models is beyond the scope of this work. 

For general beam polarisation, a nonzero mean value of a CP odd observable can 
be generated by genuine CP violation in the reaction, or by the CP odd part of the spin 
density matrix in the initial state. According to the above estimate, the latter would 
require nonstandard physics to be experimentally visible, especially for longitudinal 
beam polarisation, and thus be interesting in its own right. Similarly, a nonzero result 
for a CPT odd observable can originate from absorptive parts in the process, or from 
effects of the amplitudes with zero total helicity of the initial beams. If any such 
effects were observed, one could in a next step investigate their dynamical origin, 
using their different dependences on the beam polarisation. One possibility is to use 
that in the absence of Ar = amplitudes the cross section depends on transverse 
beam polarisation only via the product PfPf^, as seen in ([TBI) . If Pf^ = 0, the cross 
section is then independent of Pf and ip~ , whereas the interference between Ar = 
and Ar = ±1 amplitudes leads to terms with Pf cos(p~ and Pf siinp~ to the cross 
section, which can experimentally be identified via their angular dependence on (p~ . 
A possibility to search for the presence of Ar = amplitudes with only longitudinal 
beam polarisation will be discussed at the end of Sect. ^. 

Returning to purely longitudinal beam polarisation and the assumption that 
Ar = amplitudes are negligible, we remark that constraints similar to the ones 

^To be precise, one must exclude final states e'^Ve e^i'e, where nonresonant graphs contribute in 
which the initial e"*" and e~ do not annihilate. 
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Table 2: Properties of TGCs under 
parity and charge conjugation. 
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discussed above cannot be derived for C and P separately, since neutrino exchange 
maximally violates both symmetries. One can however classify the TGCs according 
to their C and P behaviour as shown in Table ^ As can be seen from the amplitudes 
in Sect. 0, terms of distinct P or C do not mix in the quadratic part of the total cross 
section 

o"2ij = / d(l)S2ij , (51) 

provided that phase space cuts are separately invariant under P and C. Under the 
same conditions the linear terms an in the integrated cross section for the couplings 
gf, g^, Rr and Xr vanish. This is due to the absence of the neutrino exchange graph 
for right handed electrons and the fact that those couplings differ in their P or C 
eigenvalues from the TGCs in the SM. Finally, real and imaginary parts of couplings 
do not mix in ij2ij- 

In the LR-basis, (J2ij has an additional block diagonality, with two separate blocks 
for the R- and the L-couplings, which cannot mix in the total cross section. Any 
block diagonality of a2ij means that already before the simultaneous diagonalisation 
(PH) the subspaces corresponding to these blocks are perpendicular to each other 
with respect to the scalar product (a, b) = J^ij o^i ^2ij bj. As a consequence, two row 
vectors of A or two column vectors of which correspond to different blocks are 
perpendicular with respect to the standard scalar product. In the case of A~'^ this 
follows from the first equation of (|4^ ) by solving for (J2, and in the case of A by 
solving for o"^^. The comparison of the matrix products AA'^ and (A~^)'^ A~^ with 
this expected orthogonality thus provides a good way to test the numerical results for 
A and A^^. Note however that the transformation described in Sect. |3.1| cannot be 
carried out on smaller blocks than those given by the four classes (a) to (d) because 
the left handed couplings mix with the right handed ones in V{0), c and V{h). 



4 Polarisation 

In this section we explore how the sensitivity of our process to anomalous TGCs 
depends on the longitudinal polarisation of the initial beams. To this end we will 
introduce an appropriate polarisation parameter P and analyse how the eigensystem 



determined by (^41) depends on it. This may be seen as a preparation for interpreting 
the numerical results in Sect. ^, which will be given in terms of the same parameter. 
To start with, we introduce a convenient notation to make explicit the polarisation 
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dependence of the matrices (3"2 and V, which have to be diagonahsed simuhaneously 
according to (@). In the following we restrict ourselves to the case of longitudinally 
polarised beams. Since in the limit rrie only two beam helicity combinations 
contribute to the amplitudes, we can use (|l^) to write the differential cross section 
S = da\p/d(j) in terms of the cross sections and for purely left and right handed 
e~ beams (and opposite e"*" helicities), 

S = P^S^ + P^S^, (52) 

where 

P^ = 1(1-P-)(1 + P+), = 1(1 + P-)(1-P+), (53) 

= (icr(_+)(_+)/(i0, = (icr(+_)(+_)/(i0 (54) 

and d(f) = d{cosQ) d{cos^) dip dicosiD) dip. Since we only deal with longitudinal pola- 
risation in this section we drop the subscript '/' in the polarisation parameters. Note 
that < P^,P^ < 1. Integrating over d(j) one obtains the total cross section as 

a = P^a^ + P V. (55) 

In the LR-basis we have 

i id 



a" 



We denote again the real or imaginary parts of the anomalous TGCs by but in 
contrast to (^) we now explicitly write an index a = L,R, so that i and j only run 
from 1 to 14. Using vector and matrix notation we can rewrite the total cross section 
as 

a = ao{l + h^a-i + hFdr2h), (58) 

where 

(To = 

^1 = — ( nR_fl, 1 > ^ 





with vectors {h°')i = and (cr"). = a^^ and the matrix (a2)ij = <J2ij = (J2ij- In the 
LR-basis and for longitudinal polarisation we thus obtain the following expression for 
the covariance matrix (IST 



ya, ^ pap, [ ^ j ^^^1 - ) , (60) 
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with a,b = L, R. Let us now investigate in detail how the eigensystem of (3"2 and V 
depends on and P^. It is useful to express P^ and P^ by new variables r and P 
(to be specified below), so that o"2, V and hence also their eigensystem only depend 
on P. For this purpose we introduce P^ and P^ through 



for some well-behaved rescaling function r{P^, P^), and define the "r-normalised" 
quantities 



0", 



' - P'^a^ + P'^a^, (62) 
SI = P^5o^ + P%«. (63) 

We then have from (|59D and (|60|): 



= ^)[ P-(P)^J' 



ya, ^ P\P)P\P){-^Jd<P^-y-^^\, (65) 



ca qb a b 

^a5(P) 7 55(P) [a5(P)]^ 

where we have made the dependence on P explicit. Note that the left hand sides 
of (^4|) and (|65D depend on P but not on r, since (T2 and do not change when 5* is 
multiplied by a constant. For the matrix M in (|i6| ) we get 



M=^^ r^_ , (66) 





where (T2 = MM is the Cholesky decomposition of the P independent submatrix 
of (J2. Then the result for the transformation (^7|) is 

X-\P)^1= ( ^_ 1 (67) 

The factors \/ P^ and \/ P^ in the rightmost matrix determine the mutual normalisa- 
tion of the blocks of left and right handed couplings. They let become singular 
in the limits P^ or P^ 0. This is not surprising because with beams of purely 
longitudinal polarisation one is sensitive to only half of the couplings. The coefficient 
((Tq(P))~^/^ in (^) leads to an overall normalisation which strongly depends on the 
polarisation. At y/s = 500 GeV we have for instance 



tLR = Ja^/a^ - 17, aS^ao^(P^ + P^/172), (6^ 
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whereas at ^/s = 3 TeV the ratio tin is about 30. From ( ^0]) we see that the matrix 
is orthogonal for any P . In the case of pure polarisation it is block diagonal 
in the left and right handed couplings. This is however not the case for general 
(longitudinal) polarisation since the diagonalisation cannot be reduced to smaller 



blocks than those given by the four discrete symmetry classes introduced in Sect. 3^ 
We now specify the transformation (|6T|) by choosing 

(69) 




and defining a polarisation parameter 




with values — 1 < -P < +1. We then have 

P^'^ = {l±Pf. (71) 

In terms of the individual beam polarisations P~ and P~^ the parameters r and P 
are given as 



1 



1 - p+p- + V (1 - P^P'Y - (P^ - P'Y j , (72) 
p - p+ 

(73) 



1 - p+p- + v(i - P^P'f - {P^ - P'f 

The reason for this particular choice is as follows. For electron polarisation Pq~ and 
positron polarisation P^ = — Pq" simply has P = Pq - Fo^' general polarisations 
P is between P~ and — P^, and the differential cross section 5* for {P~,P~^) equals 
the one for (Po~=P, P(l=—P) up to a constant. The eigenvalues c[ (cf. (^)) are 
hence the same for (P^,P+) and for (Pq",Pq^). 

To develop some intuition of how the generalised eigenvalues of (|6^ and (|65|) de- 
pend on P, we consider the case of only one left and one right handed coupling. More- 
over, we neglect the second term in (|65|), which appears only in symmetry class (a). 
The matrix W2 in (|59| ) then reduces to a single number {W2)ii, the vector (6*1)° has 
only one component s"" = {Si)i, and the 2x2 matrices which have to be diagonalised 
according to (0) can be written as 

a - ^^^^11 ( P"^ \ 

- ^)[ P« j ' ^^^^ 

1 f P^ \ f v^^ v^^ \ f P^ \ 
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.LL I Jj, \^ I „.LR _ J A ..RR ' ^ 



where 

^"'"-J'^sm- 

As in (^) we construct a symmetric matrix 

C = 7^\ i^\\ \R rrW 1' (77) 

whose usual eigenvalues are equal to the generalised eigenvalues of V, i.e. to the 
diagonal entries of the transformed matrix V in (H). They are given by 





C± = [P^V^^ + P^V^^ ± sj (P^vLL - pRyRRy + 4pLpR (^yLRy ^ _ (73) 

We approximate the matrix entries (|76D by 



v^\P)= A. / A. p > a,b = L,R, (79) 



~ab 



V 



with constants v'^^ , which should take into account their P dependence sufficiently 
well for a qualitative model. In Fig. ^ we plot the eigenvalues c± in arbitrary units 
for ctq = 1, and = (Jo/tj^Fi = (1/17)^, with the last number taken from (^). 
The ratios y^^/y^^ ^ yLRj-RR ^ jj-^ Q^j, choice of parameters are motivated by 
the power of in (^), which corresponds to the power of the neutrino exchange 
amplitude in the cross section. One can show analytically that the slopes of the 
curves for c±{P) tend to zero for P ±1. For P = 1 this cannot be seen in the 
plot, since for large Ilr (i.e. ctq S> a^) the eigenvalues change rapidly as the e~ beam 
becomes purely right-handed. To see the horizontal tangent we plot a second example 
in Fig. 1^ with a more moderate value of Ilr. Notice that for nonzero v^^ the two 
curves for c+ and c__ do not touch. If P^ or P^ is zero the matrix a"2 in ( [7^ ) and 
hence C in ( |77D is singular, which leads to a zero eigenvalue c_ at P = ±1. 
As in Sect. let 

-) (80) 

be the matrix whose columns are the normalised eigenvectors of C. We can see 
from Figs. ^ and ^ that for P = — 1 the vector £C_|_ with large eigenvalue has only an 
upper component (corresponding to h^), whereas the vector x_ with zero eigenvalue 
has only a lower component (corresponding to h^). For P = +1 the situation is 
reversed. This reflects the fact that one is only sensitive to the left handed couplings 
for P = — 1 and to the right handed ones for P = +1. We finally plot the elements 
of the transformation matrix in (^) using the notation 



ail 0,12 

0-21 0-22 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 
P P P 



Figure 6: Eigenvalues (^) and entries of the matrices (|80|) and (|8T|) for cIq = 1, 
= (1/17)2, = 1, v^^ = 0.01, 5-^^ = 0.0025 and (^2)11 = 1. 




Figure 7: Same as Fig. || but using ctq = 1, = 1/10, = 1, v^^ = 0.15, 
^flij ^ Q i g^^^ {'^2)11 = 1- 
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Writing the transformed couplings as 



we can see that for P = —1 the right handed contributions to both h'__ and h'^ vanish, 
whereas for P = +1 the same happens to the left handed ones. This behaviour of 
has to be taken into account when carrying out the simultaneous diagonalisation 
for high degrees of longitudinal polarisation, because it leads to a singularity of its 
inverse in the limit P ^ ±1. 

We have seen that under the condition that only beam helicity combinations 
with At = ±1 contribute to the cross section, the expectation values of normalised 
observables only depend on the polarisation parameter P. Such a statement no longer 
holds if beams with hehcities coupled to At = contribute as well. This provides a 
possibility to disentangle effects from genuine CP violation or absorptive parts from 
those due to nonzero At = amplitudes. In particular, the CP odd part of the spin 
density matrix for longitudinally polarised e~^e~ beams is proportional to P"*" + P~ 
and will give different contributions to normalised observables if P"*" + P~ is varied 
for fixed P. 



5 Hardly Measurable Couplings 

The particular form of the SM amplitudes for the process e~e~^ W~W^ has con- 
sequences on its sensitivity to the couplings in the CP conserving sector, which we 
shall now discuss. To this end we write the TGC part of the transition operator as 

^TGC^ E iliK. + Htm. (83) 

a=L,R i=l 

where for simplicity of notation we label the respective couplings gi, k, A, g^, g^, 
k, A by an index i = 1,...,7. Here h^^ are the SM couplings in the LR-basis 
and iff = Reiff + i Imiff are the complex anomalous couplings (which we write in 
uppercase to distinguish them from the real- valued parameters h'}). The nonzero SM 
couplings are (see (^I]) , (|l^) 

^0,1 = ^0,2 = 2^ + 4 sin^ Owil - , ^Ji = <2 = 4 sin^ 9^(1 - . (84) 

We first consider longitudinal polarisation where the differential cross section S 
(cf. |2^)) is given by (^2|). For there is no neutrino exchange contribution, so that 
we have from (|), (|^) and ( pB]) 

s"" oc E(/l^fl+-)(/l^fl+-)* 

X (/i J, + ReH^ + i ImH^) (/i£. + ReHf - i ImHf) . (85) 
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Consider now the following direction in the space of right handed anomalous cou- 
plings: 



Im H' 



where we assume ^ 1 and use the vector notation 





( 




( 


\ 






- 














1 



With (p6D the second line of (|85D equals 



R 



"-OA "-OJ 



In the space of the imaginary parts of right handed couplings there is hence a direction 
in which the differential cross section for unpolarised or longitudinally polarised beams 
has no linear term in u, but is only sensitive to order uj"^. This direction is determined 
by the real SM couplings as given by (0). Therefore, one of the functions and 
the corresponding observable O- in ( ^OD are identically zero, and V{0) contains one 
(usual as well as generalised) eigenvalue in symmetry class (b) that is zero for all 
values of P. This is confirmed by our numerical results shown in Figs. 0, |l^ and 



In the tables in Sect. |^ the eigenvalues of symmetry class (b), Cg, . . . , c'^g, are given in 
decreasing order. Using this notation we have 

S[,iM = 0, c[, = 0. (89) 



From the total rate one can derive constraints on this coupling as explained in Sect 
Now consider 







(90) 



which merely "stretches" the right handed SM couplings by a factor (1 +uj). Then 
the last line of (85) becomes 

{i+urh^hi. (91) 

In the case of purely right handed electrons or left handed positrons, i.e. for = 0, 
we have S oc from ([52|), so that the anomalous coupling (|90|) only increases the 
total rate but does not affect the normalised distribution a^^S. This holds both to 
order u and to order u^. Symmetry class (a) therefore contains a fifth zero eigenvalue 
for P = +1, in addition to the four zero eigenvalues from the left handed couplings 
Kegi, KeKi, ReA^ and R.eg^, which cannot be measured for P^ = 0. This is again 
confirmed by our numerical results (cf. Figs. |^, |13|, [l^). For P^ ^ 0, however, 5''^ also 
contributes to S. Since the functions o'^S'^ are not identical for a = L and a = R, the 
enhancement of due to ( pi]) will not just change S by an overall factor, but also 
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modify the normalised distribution cx^^S. The latter is sensitive to the anomalous 
TGC in ( pO|) in the linear approximation, since (|9lD contains a term linear in u. In 
contrast to symmetry (b), there is thus no eigenvalue which is identical to zero for all 
values of P. Note that contains interference terms of the left handed anomalous 
amplitudes and the SM neutrino exchange, so that the arguments above do not apply 
to the subspace of the left handed couplings. Also for the symmetries (c) and (d) 
there is no similar argument because CP violating TGCs are absent in the SM at tree 
level. 

The direction (^) in coupling space becomes measurable in the linear approxi- 
mation with e~^e~ beams of transverse polarisation. In fact, abbreviating 



•^0 



AT 



=+1 



{m\+-), 
{f\%''\+-), 



A7 



=-1 



(/l^o|-+), 



(92) 



where Tq is the transition operator in the SM at tree level and T"" is defined in (pSf), we 
can write the part of the differential cross section ([16D that is linear in the anomalous 
TGCs as 



(93) 



Mp ^ P'^ Re {Al*Al) ReHt - Im {Al*Al) ImH^ 



cos(^ + 2$) [Re (Ao"*AI) ReH^ - Im ImH^ 



- r sin(^ + 2$) Im Uq "^H Rei/f + Re Uq Imi/f 



where of course a = R implies r = 1 and a = L implies r = — 1. As in Sect. |^ we 
denote the degrees longitudinal and transverse polarisation of the e"^ and e~ beams 
by Pf^ and P^, and as in Sect. | we abbreviate P^ = (1 - P;")(l + P;+)/4 and 
P« = (1 + Pf)(l - P,+)/4. Note that cicr|Jj'' depends on the W' azimuthal angle $ 
only via the explicit trigonometric functions in (|93|). One thus only has da\^p'^ = if 
the three lines of (^3]) vanish separately. The first line is however the same as what we 
had for purely longitudinal polarisation, so that it vanishes for generic polarisations 
P"^ only if condition (|86| ) is fulfilled and = 0. Then relation (p3D becomes: 



da\';^ OC ooY.hU-P^'^^i'^o'^Ar) (94) 

i 

^* ^ cos(^ + 2$) Im (^Ao^*At^) + sm{^ + 2$) Re (A'* A+^J 



u ■ 



PfPt' 



cos(7/' + 2$) Im 1*^+1) + sm{ilj + 2$) Re [Ao^*A, 



where for the equality we have used the fact that Aq^ = J^ih^iAf^. Since Aq^ 
contains the neutrino exchange graph, da\^"^ no longer vanishes. Transverse beam 
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Figure 8: Schematic view of the constraints in the case of symmetry class (b). 

polarisation thus allows for the measurement of the anomalous coupling (^61), which 
is hardly possible using only longitudinal polarisation. 

Before presenting our numerical results for unpolarised and longitudinally po- 
larised beams, we must explain how to take into account the zero eigenvector (^) of 
V{0) in the analysis. From ( p3D and (|40|) we obtain for the inverse covariance matrix 
of the couplings h 

V{h)-^ = N{A-Yc'A-\ (95) 

where the hi are again 28 real parameters, viz. the real and imaginary parts of g^, 
g^, Kji, kl, etc. We number the couplings in the order of their symmetry class (a) 
to (d), and within each symmetry class take the L-couplings first. We then have 
hi3 = Img^ and = Im/t/j. Note that V{h)~^ always exists, even in our case 
where one parameter is unmeasurable. In this case V{h)~^ is a singular matrix with 
a one- dimensional zero eigenspace coming from c[q = 0. Geometrically speaking, 
the error ellipsoid defined by V{h)~^ is degenerate in such a way that the length 
of one principal axis is infinite. Instead of an ellipsoid we have a cylinder whose 
axis corresponds to the direction of the unmeasurable coupling and whose cross- 
section (orthogonal to the axis) is an ellipsoid giving the errors on the couplings in 
the orthogonal space. We know from ( plD and (^) that the unmeasurable direction 
is given by Img^ = Iirkr ^ with all other couplings being zero. Therefore the 
projection of the cylinder onto the (Im^ff )-(Im Kj:j)-plane is a band in the lm{gi+K,{i)- 
direction, see Fig. |^(a). This shows that we cannot obtain any constraint on Img^ or 
Im unless one of them is known. We can however choose coordinate axes parallel 
and orthogonal to the band shown in Fig. H(a). In other words, we perform a rotation 
by —45° in the (Imgff )-(Im Kij)-plane, 

h = Rh, (96) 
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where R is the identity matrix except for the {Im gi)-{lm K,pi)-hlock, which reads 

The new couphngs hi are the same as the hi, except for hi^ = Im((7f — kr)/\/2 
and flu = Im(5ff + Kfi)/\^, which replace Img^ and ImKR. The inverse covariance 
matrix of the new couphngs is 

vCh)'^ = RV{h)-^R^. (98) 

Ah entries in the 14th row and in the 14th column of V{h)^^ are equal to zero: there is 
no correlation between the unmeasurable Im^g^+Kn) / A/2-direction and the couplings 
hi with i 7^ 14. These couplings are hence constrained by a 27-dimensional ellipsoid, 
which is drawn schematically in Fig. ||(b) for one further coupling, taken to be Im A/j. 
This ellipsoid is determined by the "reduced" 27 x 27 matrix V^~}{h) obtained from 
V{h)~^ by deleting the 14th row and the 14th column. Its inverse Vred{h) is the 
covariance matrix of lm{g^ — kr) / and the other 26 measurable couplings. In 
particular, the width of the band in Fig.JKa) gives the error on Im((7f — kr)/\/2 
in the presence of all other 27 couplings hi, cf. Fig. ^(b). We finally mention that 
because of the discrete symmetries explained in Sect. p.3| , the matrix V{h)~^ is block 
diagonal with one block for each symmetry class (a) to (d), so that the errors on the 
couplings of class (a), (c) and (d) are entirely unaffected by the previous discussion. 



6 Numerical Results 

In this section we present the results for the generalised eigenvalues c- of the covari- 
ance matrix V{0) and the corresponding errors 5h[ = {Nc'i)~^^'^ on the transformed 
couplings. The covariance matrix for the couplings in any other parameterisation 
is then obtained by conventional error propagation. We discuss its most important 
properties in the LR-basis for -y/i = 500 GeV and unpolarised beams in Sect. |0| . In 
Sect. |6.2| we investigate the gain in sensitivity by longitudinal e~ as well as additional 
polarisation. The results for higher cm. energies are reported in Sect. p.'S\ In 
Sect. ^]4| we finally give the constraints which can be obtained from the total rate 
according to (|4l|) . Numerical rounding errors on the results presented in this section 
are typically of order 1%. 



We use myy = 80.42 GeV and mz = 91.19 GeV from ||25[, and the definition 
siia^ 9\Y = 1 — m^/m^ ^r the weak mixing angle. For the total event rate N of 
the semileptonic channels with e and /i summed over we use the values listed in 
Table ^. They correspond to an effective electromagnetic coupling constant a = 
1/128 and integrated luminosities of 500 fb~^, 1 ab~^ and 3 ab^^ at y/s = 500 GeV, 
800 GeV and 3 TeV, respectively. We assume full kinematical reconstruction of the 
final state, except that the jet charges are not tagged. Due to this twofold ambiguity 
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Table 3: Total event rate N in units of 10^ for the semileptonic channels with e 
and II summed over. Corresponding luminosities are given in the text. P~ and 
respectively denote the degrees of longitudinal polarisation of the e~ and e"*" beams, 
and P is given in ([73|). 



polarisation 


[GeV] 


P- 


P+ 


P 


500 800 3000 


-80% 


+60% 


-71% 


3280 3410 1280 


-80% 





-50% 


2050 2130 799 











1140 1190 444 


+80% 





+50% 


235 242 89.7 


+80% 


-60% 


+71% 


103 103 37 



we cannot take the zeroth and first order parts of the differential cross section 5* 
as the denominator and the numerator of the optimal observables O, but use their 
respective sums over the two experimentally undistinguished final states (cf. [0). 

6.1 Unpolarised Beams at 500 GeV 

We first consider the sensitivity at ^/s = 500 GeV with unpolarised beams. In Ta- 
bles we list the standard-deviation 5hi = [V{h)iiY^'^ for each coupling hi, which 
gives its error in the presence of all other couplings. Notice the difference of this to 
(A^Cjj)"^/^, which corresponds to the error on hi when all other couplings are assumed 
to be zero. We also give the correlation matrix 

W{h),, = . ^^^^'^ (99) 
^Vih)uV{h),j 

of the couplings for each symmetry class (a) to (d). In the case of symmetry (b) we 
use the reduced matrix Ked(^) introduced in Sect. ^. Since W{h) is symmetric we 
only list its upper triangular part. 

The values of 6h range from about 5 ■ 10~^ to 10~^ within each symmetry class. 
The smallest are those for X^, ni and since at high energies the corresponding 
terms in the helicity amplitudes contain a factor 27^ (cf. Table |l]). In all cases the 
errors on the R-couplings are larger than those on the respective L-couplings, viz. by 
a factor 1.5 to 3.4 for Im k, Re Ak and Re g^, and by a factor between 4 and 7 for the 
other couplings. This is because (for unpolarised or longitudinally polarised beams) 
the //-exchange interferes with the amplitudes containing the h^, but not with those 
containing the h^. In general, the sensitivity to the real part of a specific coupling 
is roughly of the same size as the sensitivity to its imaginary part, the errors on the 
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latter being rather larger for CP conserving couplings and smaller for the CP violating 
ones. To get an accurate picture of the sensitivities, correlations have to be taken into 
account. Looking at the 2x2 blocks corresponding to hf and for a given index i, 
i.e. at the diagonal entries of the right upper block in the correlation matrices, we 
see that the absolute values of the correlations are smaller than about 0.2, except 
for ReA/t, Imn and Img^, where they are still smaller than 0.6. The corresponding 
correlations would be substantial in the basis of the 7- and Z-iorm factors, which is 
hence not a very suitable parameterisation for the present reaction, compared with 
the LR-basis (cf. [0). Considering the matrix blocks of correlations among different 
L-couplings or among different R-couplings, we find that about half of them have an 
absolute value larger than 0.4. Note that there are correlations of order 0.5 between 
couplings with different C or P eigenvalues. 

By simultaneous diagonalisation (see Sect. p.l|) we determine the generalised eigen- 
values and corresponding errors 6h'^ given in Tables ^ and ^. For symmetry class (b) 
we give the transformation matrix A and its inverse in Tables 10 and In our 



numerical calculation we find the smallest eigenvalue in symmetry class (b) to be 
c[q ~ 10"^'*. We can however use the result of our analytical considerations in Sect. ^ 
and set c[g to zero. The same holds for the last column of A except for its Img^- 
and Im Kij-components, which determine the corresponding "blind" direction in the 



LR-basis. As explained in Sect. |3.3| the P odd coupling (75 does not mix with the other 
couplings in (J2, and the same is true for the left and right handed couplings. From 
this block structure of (72 and the relation = {A~^)^ it follows that in the last 
row of A^^ we can set all entries to zero, except for the Img^-, ImKji- and ImA/j- 
components. Numerically we find that the absolute values of those matrix entries 
which we set to zero are smaller than 10~^ for A and smaller than 10~^ for A~^. We 
remark that we have computed the matrix A~^ by inverting A using singular value 
decomposition As mentioned at the end of Sect. we have as a further check 
evaluated the products AA^ and (A^^)'^ A'~^ . 

6.2 Polarised Beams 

At future e"'"e~-colliders longitudinal polarisation of both initial beams is envis- 



aged |3^, An electron polarisation of P = ±80% and a positron polarisation 



of P"*" = ±60% is considered to be achievable. 



In Tables |T2| and |T3| we give the errors 6h on the real couplings (in the presence 
of all couplings) for ^/s = 500 GeV and various combinations of beam polarisations. 
For all couplings and all couplings we find roughly the following gain or loss 
in sensitivity using always the event rates of Table ^. Turning on e~ polarisation of 
—80% we gain a factor of 1.4 for and loose a factor of 6 for h^. If in addition 

^The matrices A^^ for symmetry classes (a), (c) and (d) can be obtained from the authors. 
Further numerical results for unpolarised and longitudinally polarised beams at various cm. energies 
are given in BQ]. 



31 



= +60% we gain a factor of 1.8 for and loose a factor of 17 for compared to 
unpolarised beams. For P~ = +80% we loose a factor of 2.6 for and gain a factor 
of 3.0 for h^. If furthermore P+ = —60% we loose a factor of 5 for and gain a 
factor of 5.5 for compared to unpolarised beams. Especially for the right handed 
couplings the gain from having both beams polarised is thus appreciable. 

The behaviour of the generalised eigenvalues as a function of the parameter P 
introduced in Sect. ^ is shown in Figs. p|-pO[ Although the four largest eigenvalues 



are more or less constant for P < 0, the transformation matrix A ^ is not. This can 



be seen from Tables |T^, ^ and |T6|. For the largest eigenvalue c[ of symmetry (a) we 
find that the smaller P is, the more the R-components are suppressed, i.e. the more 
one purely measures the h^. Going from P = to P = 1 we become more and more 
sensitive to the h^. For the fourth lowest curve in Fig. |^, corresponding to Cg, as well 
as for the smallest eigenvalue Cg of symmetry (a) we find the opposite tendency. Note 
that in the case of ±100% electron or positron polarisation we can only be sensitive to 



at most half of the couplings. This is seen for symmetries (c) and (d) in Figs. |TT|, 12 
15] , |T6|, |19| and |2y: half of the curves go to zero at P = ±1. For class (a) (cf. Figs. 



|T3| and we find one additional eigenvalue going to zero at P = +1 and for class (b) 
(cf. Figs. |T1| and 0) there is a zero eigenvalue for all P, as explained in Sect. ^. 
Comparing with Fig. ^ we see that for symmetries (b) to (d) the shape of the curves 
is qualitatively well described by the simple model of Sect. Although the lower 
and upper curves for c± do not intersect in our examples there, an intersection like in 
Fig. |ll] is not excluded. In general, it is however not possible to associate a certain 
pair of couplings to a pair of curves in Figs. p|-pO| for the full range of P. This is 



particularly obvious from the eigenvalues of symmetry (c) at = 3 TeV (Fig. ^), 
where some curves alternately play the role of the lower-type and upper-type curves 
in the simplified model. Moreover, for symmetry class (a) the description of the shape 
of the eigenvalue curves is less obvious due to the second term in the brackets of (^) . 

6.3 Energy Dependence 

The gain in sensitivity when going up from 500 GeV to 800 GeV — using always the 
event rates of Table ^ — lies between 1.4 and 2.7 for all couplings except for Im^^, 
where it is 3.6. At 3 TeV we gain a factor of about 25 compared to 800 GeV for this 
coupling, and of 1.5 to 8 for all others. For symmetries (a) and (c) we give 6h in 



Tables 17 and 18 



Note that this gain is not due to the the total rate, which actually decreases with 
energy (cf. Table |^). The largest gains are achieved for k, A and A, which have a 
prefactor 27^ in the amplitude. We remark that both for real and imaginary parts 
the gains in sensitivity for an L-coupling and the corresponding R-coupling are of the 
same size, except for ImKR. Furthermore, except for Akl, g^, gf and Rr, the gain 
is slightly larger for the imaginary than for the real parts. For the real parts of the 



couplings we also give the errors on the transformed couplings 5h'^ in Tables |l^ and ^ 
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Note that the transformations (|37|) are not identical at the various cm. energies, and 
neither are the couphngs h[. Due to the different normahsation of the h'^ achieved 
by (pTl) their errors dh'^ may well increase with rising energy although the errors in a 



fixed basis as in Tables 17 and 18 decrease. 



6.4 Constraints from the Total Rate 

As explained in Sect, ^the measurement of the total cross section restricts the anoma- 
lous TGCs in the /i^-basis to a shell between two hyperspheres in the multi-dimensional 
parameter space. For the couplings given in a basis before the transformation we have 
hyperellipsoids instead of hyperspheres. With y/s = 500 GeV and unpolarised beams 
the expansion of the total cross section (|41| ) is 

a/ao = 1 -0.026 (100) 



+ {h[ + 0.16)2 + (/i2 + 0.026)2 + (/ig + 0.0042)2 ^ ^^^^ ^ 0.0061)^ 
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+ (/i5 - 0.013)2 + [Hq - 0.022)2 ^ ^^/^ _ 0.0093)2 ^ ^^^/^ ^ 0.00013)2 + ^(/i-)2. 

i=9 

At 800 GeV we obtain 

a/ao = 1 - 0.016 (101) 
+ {h[ + 0.13)2 + {h'^ + 0.0078)2 ^ ^^/^ ^ 0.0025)2 + {h'^ + 0.0027)2 

28 

+ (/I's - 0.0066)2 ^ ^j^f^ _ 0.013)2 + {h'^ - 0.0062)2 ^ ^j^f^ _ 0.00023)2 + ^(/i-)2, 

i=9 

and, finally, at 3 TeV 

a/ao = 1 - 0.0071 (102) 
+ {h[ + 0.084)2 ^ (^^/^ ^ 0.00083)2 ^ ^^^/^ ^ 4 3 . ^g-6^)2 ^ ^^^^ ^ 0.00060)2 

28 

+ (/I's + 0.0030)2 _^ (^^/^ _ 0.0052)2 ^ ^^/^ _ o.ooi7)2 + {h'^ + 3.2 ■ 10^^)2 + ^(/i-)2. 

1=9 

We remark again that the couplings /i^ are not the same at different energies. For 
a measurement of the rate with a (purely statistical) error ^/N the thickness of 
the shell is 5.8 ■ 10"=^ at 500 GeV, 7.2 ■ 10"^ at 800 GeV and 18 ■ lO'^ at 3 TeV 
(cf. Table |^). Systematic errors could be more important. The results ( [100|) agree 
quite well with those of |l^ for all couplings except for the smallest term h'g and 



for /13. Note that these results strongly depend on a reliable transformation matrix 
A. In numerical instabilities occurred in the diagonalisation procedure, whereas 
here A is obtained iteratively as explained in Sect, ^and was found to be stable. 
It has been pointed out ]18[ that the constraints from the total rate are in general 



of the same size as the largest error on the couplings determined from the normalised 
distribution, which we confirm. 
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7 Conclusions 



In this paper we have investigated the sensitivity of optimal observables to anomalous 
TGCs in e~e+ W~W^ /1/2/3/4 at future linear colliders. We have treated all 
28 couplings at a time and disentangled them by a simultaneous diagonalisation of 
the covariance matrix of the observables and the term of the integrated cross section 
which is quadratic in the anomalous couplings. Several relations arising from discrete 
symmetry properties of the differential and the total cross section, based on the 
classification of the TGCs according to their CP and CPT eigenvalues, have been used 
to simplify the analysis and to check the numerical stability of the results. We have 
re-examined the conditions for investigating CP violation in the presence of polarised 
e~^e~ beams. Suitably defined CP odd observables receive contributions only from 
genuine CP violation in the interaction and from the CP odd part of the spin density 
matrix of the initial e+e~ state. The two kinds of effects can be separated using their 
different dependence on the beam polarisations, either with transversely polarised 
beams (cf. Sect. |3.3| ) or with suitable combinations of longitudinal polarisations (cf. 
Sect. ^). An analogous statement holds for the investigation of absorptive parts in 
the scattering amplitude by CPT odd observables. 

We find that already at a 500 GeV collider with unpolarised beams and the event 
rates of Table 0, the statistical errors obtained with the optimal-observable method, 
treating all couplings at a time, are considerably smaller than the constraints ob- 
tained from single parameter fits of LEP2 data. Our errors for ^/s = 500 GeV and 
800 GeV treating all couplings are of the same size or smaller than those obtained at 
generator level for TESLA by a spin density matrix method with a restricted number 
of couplings 0, |32| . 

We have performed a detailed study of the sensitivity to anomalous TGCs for 
different longitudinal polarisations and cm. energies of the lepton beams. The pola- 
risation dependence is conveniently expressed through a parameter P. A simple 
model has been analysed (see Figs, ^, which provides an understanding of the 
dependence of sensitivities on P. We find that beam polarisation can provide an 
appreciable gain in sensitivity, especially for right handed couplings (cf. Tables [T2| 
and |l^). At y/s = 500 GeV the sensitivity to these couplings increases by a factor 
of about 3 when going from unpolarised beams to -|-80% electron polarisation. With 
additional —60% positron polarisation this factor is about 5 or larger. 

In the sector of CP conserving imaginary TGCs we have found one linear combi- 
nation of couplings which appears only quadratically in the differential cross section 
for longitudinally polarised beams (cf. Sect. |^). Therefore the normalised event dis- 
tribution of W pair production for unpolarised or longitudinally polarised beams does 
not provide a good way to measure this coupling — regardless of whether the analy- 
sis uses optimal observables or any other method. Information on this coupling can 
however be obtained either from the total event rate, or from the normalised event 
distribution with transversely polarised beams. 
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Since our numerical results are at the "theory level", a study using Monte Carlo 
event samples and including a detector simulation will at some point be necessary. 
As noted in Sect. ^ inclusion of detector resolution and phase space cuts is easy in 
the framework of the method presented here. 

In this work we have analysed the normalised event distributions in the linear 
approximation for anomalous TGCs. We emphasise that the extension of our method 
to the fully non- linear case is available [|18[. Such an analysis for the non- linear case 
has indeed been performed successfully with LEP2 data . 

To summarise, we have performed a detailed study of the sensitivities achievable 
in the measurements of anomalous 'yW~^W~ and ZW^W~ couplings at future linear 
e+e" colliders. We advocate the use of integrated optimal observables for this purpose. 
We have shown that our method allows for a very good overview of the sensitivities 
and of their dependence on the beam polarisations in the space of the 28 anomalous 
couplings. Such measurements will be a crucial check wether the structure imposed 
by the fundamental gauge group SU(2)l x U{1) of the electroweak interactions in the 
Standard Model for the couplings of three gauge bosons is realised in nature. 
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Table 4: Errors 5h in the presence of all other couplings and correlation matrix W{h) 
at ^/s = 500 GeV with unpolarised beams for the couplings of symmetry (a) (see 
Sect. p.3D , i.e. for the real parts of the CP even couplings. 



h 




ReAgf 


Re Akl 


ReAi 


Reffi^ 


Re Agf 


Re Akii 


Re Afl 


Re 


ReAgf 


2.6 


1 


-0.60 


-0.35 


0.21 


-0.070 


0.25 


-0.054 


0.15 


ReAuL 


0.85 




1 


0.096 


-0.054 


0.20 


-0.59 


0.13 


0.019 


ReAi 


0.59 






1 


-0.034 


0.099 


-0.080 


0.030 


0.10 


Re 


2.0 








1 


-0.084 


0.11 


-0.13 


0.075 


ReA^f 
Re Akr 
ReXji 
Regi 


10 
2.4 
3.6 
6.7 




1 


-0.70 
1 


-0.56 
0.023 
1 


0.65 
-0.34 
-0.25 
1 



Table 5: Same as Table ^ but for symmetry (b), i.e. the imaginary parts of the CP 
even couplings. As explained in Sect, ^we have 5Im((yff + = oo and no correlation 



of this coupling 
h 


with the 

Shx 10^ 


others. Thus we only give th 

Im gi Im kl Im A l Im .95 


e reduced 7x7 matrix here. 

■^lm{g^ - kr) IitiXr lmg§ 


Im Kl 
Im Al 
Imyf 


2.7 
1.7 

0.48 
2.5 


1 -0.47 -0.50 -0.12 
1 0.0070 0.41 
1 -0.15 
1 


0.028 0.16 0.038 
0.33 -0.10 0.68 
-0.00069 -0.21 -0.22 
0.081 0.22 0.50 


^Im(5f - kr) 

Im A_R 

Imgf 


11 
3.1 
17 




1 -0.53 0.60 
1 -0.11 
1 



Table 6: Same as Table ^, but for symmetry (c), i.e. the real parts of the CP violating 
couplings. 



h 


ShxlO^ 


Re 54 


ReAL 


Re Kl 


Re^f 


Re Xr 


Re Kr 


Reff4 


2.5 


1 


-0.055 


-0.49 


-0.091 


-0.18 


0.073 


ReAi 


0.60 




1 


0.27 


0.073 


0.0088 


-0.16 


Re 


2.7 






1 


0.036 


0.11 


0.14 


RegJ 
Re Afl 
Re kr 


10 
3.8 
11 




1 


-0.24 
1 


-0.47 
0.65 
1 
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Table 7: Same as Table ^ but for symmetry (d), i.e. the imaginary parts of the CP 
violating couplings. 



h 


(5/1x10^ 


Imgf 


Im Al 


Im 


Imgf 


ImAfl 


ImKfl 




1.9 


1 


-0.059 


0.092 


0.20 


0.22 


-0.017 


Im Al 


0.46 




1 


0.53 


-0.15 


-0.18 


-0.015 


Im 


2.0 






1 


-0.33 


-0.099 


0.14 


Imgf 
Im Afl 
ImKfl 


7.7 
2.9 
8.6 




1 


-0.12 
1 


-0.68 
0.56 
1 



Table 8: Generalised eigenvalues of the covariance matrix V{0) and the corre- 
sponding errors 5h^ on the transformed couplings, obtained from (^) and Table |^ at 
500 GeV with unpolarised beams for symmetries (a) and (b). 



i 




5h[ X 10^ 


i 


< 


5h[ X 10^ 


1 


1.44 


0.780 


9 


1.27 


0.831 


2 


1.17 


0.866 


10 


1.01 


0.931 


3 


0.751 


1.08 


11 


0.791 


1.05 


4 


0.557 


1.25 


12 


0.287 


1.75 


5 


0.318 


1.66 


13 


0.0584 


3.88 


6 


0.108 


2.85 


14 


0.0221 


6.30 


7 


0.0366 


4.90 


15 


0.0102 


9.29 


8 


0.0147 


7.72 


16 





oo 



Table 9: Same as Table ^ but for symmetries (c) and (d). 



i 




6h[ X 10^ 


i 


c[ 


5h[ X 10^ 


17 


1.17 


0.868 


23 


1.40 


0.792 


18 


0.585 


1.23 


24 


1.02 


0.929 


19 


0.320 


1.66 


25 


0.829 


1.03 


20 


0.0645 


3.69 


26 


0.219 


2.00 


21 


0.0262 


5.78 


27 


0.0316 


5.27 


22 


0.0131 


8.18 


28 


0.0241 


6.04 
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Table 10: Coefficient matrix A x 10^ for symmetry (b) 



h'. 


h'lo 


h'n 


h[2 


Ks 


h'u 


h[5 




140 

-4.3 


-85 
3.4 


-130 
-1.3 


97 

-64 


12 

-15 


7.1 

0.80 


0.86 

12 








Ini kl 
Ini Xl 



6.6 
120 



46 
-52 



10 
150 



-6.1 
-7.9 



0.87 
-16 



-1.9 
13 



-0.94 
12 



Imgf 

Im KR 
Im A_R 



6.3 -11 

-0.89 1.7 

-1.8 2.3 

-22 16 



8.2 
-1.4 
-1.2 

1.4 



-23 
5.1 
-0.50 
55 



170 

-37 
0.62 
-86 



-91 
7.6 
47 
24 



110 

-21 
-6.5 
170 



-50 

-50 





Table 11: Coefficient matrix A ^ x 10^ for symmetry (b) 



i 




lm.KL 


Im Al 






ImKR 




Img^ 


9 


42 


45 


110 


29 


0.49 


-0.49 


-6.7 


-5.3 


10 


1.5 


-12 


210 


-12 


-1.2 


1.2 


7.3 


3.9 


11 


-33 


-54 


-15 


37 


1.1 


-1.1 


-2.6 


0.35 


12 


-0.73 


-140 


2.0 


-1.9 


-5.0 


5.0 


-14 


13 


13 


-2.1 


-38 


7.2 


-4.0 


36 


-36 


87 


-21 


14 


1.3 


4.6 


-6.6 


3.1 


2.2 


-2.2 


210 


5.7 


15 


4.5 


34 


-1.5 


2.8 


19 


-19 


17 


41 


16 














-35 


-160 


-40 






Tabic 12: Errors 5hxlO^ on the couplings of symmetry (a) at 500 GeV for different 
initial beam polarisations. 



p- 


P+ 




Re Akl 


ReAi, 




ReAgf 


Re A/tij 


ReAij 


Re 


-80% 


+60% 


1.5 


0.47 


0.34 


1.1 


169 


40 


57 


112 


-80% 





1.9 


0.60 


0.43 


1.5 


62 


14 


21 


41 








2.6 


0.85 


0.59 


2.0 


10 


2.4 


3.6 


6.7 


+80% 





6.9 


2.3 


1.5 


5.3 


3.5 


0.83 


1.2 


2.3 


+80% 


-60% 


13 


4.5 


2.8 


10 


2.0 


0.47 


0.67 


1.3 
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Table 13: Same as Table |T2], but for symmetry (c). 



p- 


P+ 




RcAl 


Re 


Reg« 


Re Ai? 


Re Kfl' 


-80% 


+60% 


1.4 


0.34 


1.5 


174 


61 


193 


-80% 





1.8 


0.43 


1.9 


62 


22 


69 








2.5 


0.60 


2.7 


10 


3.8 


11 


+80% 





6.5 


1.5 


6.9 


3.2 


1.3 


3.7 


+80% 


-60% 


13 


2.9 


13 


1.8 


0.70 


2.0 



Table 14: Vector components belonging to the largest eigenvalue c[ of symmetry (a) 
at 500 GeV for different longitudinal polarisations (two digits precision). Each line is 
the first row of A'^ x 10^ in the LR-basis. 



p- 


P+ 


ReAg^ 


Re Akl 


ReAi 




ReA^f 


RbAkr 






-80% 


+60% 


-34 


-150 


-33 


18 


-0.19 


-1.2 


-0.089 


0.11 


-80% 





-34 


-150 


-33 


17 


-0.72 


-4.5 


-0.36 


0.41 








-32 


-150 


-29 


11 


-9.0 


-48 


-7.1 


4.8 


+80% 





-13 


-84 


-10 


-3.9 


-91 


-390 


-110 


41 


+80% 


-60% 


-5.5 


-42 


-4.9 


-2.6 


-190 


-840 


-240 


79 



Table 15: Same as Table M, but for the eigenvalue Cg of symmetry (a) 



p- 


P+ 


ReAgf 


Re Akl 


ReAi 




Re Agf 


Re Akii 


Re Xji 


Re 


-80% 


+60% 


-1.5 


-7.9 


1.1 


-2.1 


7.0 


14 


8.3 


-5.6 


-80% 





-2.9 


-16 


2.1 


-4.1 


13 


28 


16 


-10 








-9.5 


-50 


6.9 


-14 


37 


82 


44 


-27 


+80% 





-13 


-41 


150 


-20 


14 


42 


-9.6 


-39 


+80% 


-60% 


-20 


-72 


100 


-24 


26 


120 


-28 


-18 



Table 16: Same as Table 0, but for the smallest eigenvalue Cg of symmetry (a) 



p- 


P+ 


ReAgf 


ReAKL 


ReAi 


Re 


ReA^f 


Re Akj^ 


ReAfl 


Regf 


-80% 


+60% 


0.029 


1.0 


0.071 


-0.16 


0.81 


- 4.1 


38 


-0.91 


-80% 





0.054 


1.8 


0.17 


-0.31 


1.6 


-7.5 


73 


-1.7 








0.15 


4.8 


0.64 


-1.0 


5.4 


-21 


220 


-4.8 


+80% 





0.26 


14 


1.2 


-3.2 


17 


-61 


650 


-13 


+80% 


-60% 


-3.2 


47 


-5.1 


-8.5 


38 


-170 


1100 


-24 



39 



Table 17: Errors 6hxl0^ on the couplings of symmetry (a) for different cm. energies. 



[GeV] 




ReAKL 


ReAi 




ReA^f 


Re Akji 




Re 


500 


2.6 


0.85 


0.59 


2.0 


10 


2.4 


3.6 


6.7 


800 


1.6 


0.35 


0.24 


1.4 


6.2 


0.92 


1.8 


4.8 


3000 


0.93 


0.051 


0.036 


0.88 


3.1 


0.12 


0.36 


3.2 



Table 18: Same as Table ^ but for symmetry (c). 



V^s [GeV] 


Reff4 


ReAi 


Re kl 


Reg^ 


Re Ait 


Re kji 


500 


2.5 


0.60 


2.7 


10 


3.8 


11 


800 


1.7 


0.24 


1.8 


6.5 


1.8 


6.8 


3000 


0.90 


0.036 


0.97 


3.4 


0.36 


3.2 



Table 19: Errors 6h^x 10^ on the transformed couplings of symmetry (a) at different 
cm. energies. 



i 


500 GeV 


800 GeV 


3 TeV 


1 


0.780 


0.765 


1.26 


2 


0.866 


0.841 


1.35 


3 


1.08 


1.16 


2.02 


4 


1.25 


1.26 


2.39 


5 


1.66 


1.83 


4.18 


6 


2.85 


3.07 


5.29 


7 


4.90 


4.96 


8.54 


8 


7.72 


9.27 


20.8 



Table 20: Same as Table ^ but for symmetry (c) 



i 


500 GeV 


800 GeV 


3 TeV 


17 


0.868 


0.832 


1.35 


18 


1.23 


1.22 


2.03 


19 


1.66 


1.58 


2.52 


20 


3.69 


3.39 


5.12 


21 


5.78 


5.54 


8.74 


22 


8.18 


9.53 


20.9 
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P 



Figure 9: Generalised eigenvalues of the correlation matrix V{0) for the couplings 
of symmetry class (a) (cf. Sect. |3.3|) at ^/s = 500 GeV. The c[ do not depend on 
the total rate A^. Errors on the transformed couplings h[ ( ^7\) are obtained as 6h^ = 
{Nc'j)'^^"^. Vertical lines mark the five cases investigated in detail in Sect. |6.2| , i.e. from 
left to right (P-,P+) = (-80%, +60%), (-80%, 0), (0,0), (+80%, 0), (+80%, -60%). 
P is given in terms of P~ and by ([73D. 
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Figure 10: Same as Fig. ^ for symmetry class (b). 
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Figure 11: Same as Fig. ^for symmetry class (c). 
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Figure 12: Same as Fig. ^ for symmetry class (d). 
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Figure 13: Same as Fig. ^ for ^/s = 800 GeV (symmetry class (a)). 
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Figure 14: Same as Fig. ^for ^/s = 800 GeV and symmetry class (b). 
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Figure 16: Same as Fig. ^for ^/s = 800 GeV and symmetry class (d). 
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Figure 18: Same as Fig. |^ for y/s = 3 TeV and symmetry class (b). 
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Figure 20: Same as Fig. |^ for y/s = 3 TeV and symmetry class (d). 
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A Appendix: Conventions 

Momenta and helicities of incoming and outgoing particles are denoted as in Fig. |l|. 
We evaluate the production amplitude in the frame obtained from the one in Fig. |I| 
by a rotation of G around the y-axis, so that the new z'-axis points along the W~ 
momentum. For the respective polarisation vectors e;^ and ej of W~ and we 
choose in this frame 

e± = ^ (0,Tl,-i,0), 

eo = ^(g^O,0,g°), 
1 . 

e± = (0,Tl,i,0), 

eo = ^(-g^O,0,g°). 
mw 

The four-spinors for the initial leptons are expressed through two-spinors x iii the 
usual way [0], with 

/ e e\ / e e\ , , 

Xr=+i = (cos — , - sm — 1 , Xt=-i = ( sm ^> cos — 1 (103) 

for the electron and 

X-=+i = ( sm — , cos — 1 , X-=-i = ( - cos — , sm — 1 (104) 

for the positron. The evaluation of the diagrams in Fig. ^ then leads to ([T8|) , where 
the d-functions are defined in the usual fashion: 

<±i = ^(lircosG), (106) 
dl^2 = ±^(1 ±rcose)sine. (107) 
For the W decay tensors (|11|) one has 

dVy^ = 2AnmwT{W- ^ fJ'2)lx'lld{cos^)dip, 

dVj,j = 2AT^mw T{W+ fjl) Ij, Ij rf(cos?9) dTp , (108) 

where 

I. = d+{^)e-'^, L = rf+(^)e'^, 
^0 = -c^o(^), ^0 = -do{'&), 

1+ = d^{d)e''^, 1+ = d_{§)e-''^ (109) 
with d±{x) = (1 ± cosx)/\/2 and dQ{x) = sinx. 
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